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Abstract 

We study the large distance expansion of correlation functions in the free massive Majorana 
theory at finite temperature, alias the Ising field theory at zero magnetic field on a cylinder. We 
develop a method that mimics the spectral decomposition, or form factor expansion, of zero- 
temperature correlation functions, introducing the concept of "finite-temperature form factors" . 
Our techniques are different from those of previous attempts in this subject. We show that an 
appropriate analytical continuation of finite-temperature form factors gives form factors in the 
quantization scheme on the circle. We show that finite-temperature form factor expansions are 
able to reproduce expansions in form factors on the circle. We calculate finite-temperature form 
factors of non-interacting fields (fields that are local with respect to the fundamental fermion 
field). We observe that they are given by a mixing of their zero-temperature form factors and 
of those of other fields of lower scaling dimension. We then calculate finite-temperature form 
factors of order and disorder fields. For this purpose, we derive the Riemann-Hilbert problem 
that completely specifies the set of finite-temperature form factors of general twist fields (order 
and disorder fields and their descendants). This Riemann-Hilbert problem is different from the 
zero-temperature one, and so are its solutions. Our results agree with the known form factors 
on the circle of order and disorder fields. 
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1 Introduction 

Quantum field theory (QFT) at finite temperature is a subject of great interest which has 
been studied from many viewpoints (see, for instance, ^). The main goal of QFT is the 
reconstruction of a set of correlation functions in which the physical information provided by 
the theory is embedded. For instance, two-point correlation functions are related to response 
functions, which can be measured and which provide precise information about the dynamics 
of a physical system at thermodynamic equilibrium. For the purpose of comparison with 
experiment, it is important to know the influence of a non-zero temperature on correlation 
functions. In particular, both static (equal-time) and dynamical two-point correlation functions 
at finite temperature still need more accurate study. 

In recent years, thanks to advances in experimental techniques allowing the identification and 
study of quasi-one-dimensional systems (see for instance [3 E]), there has been an increased 
interest in calculating correlation functions in 1+1-dimensional integrable models of QFT (for 
applications of integrable models to condensed matter systems, see for instance the recent review 
Integrable models are of particular interest, because in many cases, matrix elements of local 
fields in eigenstates of the Hamiltonian, or form factors, can be evaluated exactly by solving an 
appropriate Riemann-Hilbert problem in the rapidity space jSl 13 El El • At zero temperature, 
correlation functions are vacuum expectation values in the Hilbert space of quantization on the 
line. A useful representation of zero-temperature correlation functions is then provided by their 
form factor expansion (or spectral decomposition): 



where r is the Euclidean time and r = \/a;^ + r^. This expansion is obtained by inserting 
the resolution of the identity in terms of a basis of common eigenstates of the momentum 
operator and of the Hamiltonian, parameterized by the rapidity variables 6'j's and by the 
particle types e^-'s. Here M(e) is the mass of a particle of type e (and we took the basis of 
m-states). Using many of such insertions, form factor expansions can be obtained for multi- 
point correlation functions as well. This is a useful representation because it is a large-distance 
expansion, which is hardly accessible by perturbation theory, and which is often the region of 
interest in condensed matter applications. Also, form factor expansions in integrable models at 
zero temperature have proven to provide a good numerical accuracy for evaluating correlation 
functions in a wide range of energies, and combined to conformal perturbation theory give 
correlation functions at all energy scales (an early work on this is jlUj). 

In this paper we study the large-distance expansion of correlation functions in the free massive 
Majorana field theory on a line at finite temperature. This model of 1+1-dimensional quantum 
field theory occurs, for instance, as the off-critical scaling limit of the finite-temperature Ising 
quantum chain (with an infinite number of sites) in a transverse magnetic field, near to its 
quantum critical point. It also occurs as the off-critical scaling limit of the statistical Ising 
model on a cylinder, near to its thermal critical point; the circumference of the cylinder is 
the inverse temperature of the Majorana theory. For an exposition on these scaling limits and 



(vac|C»(x,r)O(0,0)|vac) = ^ ^ 



g-r-E,M(e,)cosh(e,) 



fc=0 ei, 
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for references of early works, see, for instance, the book ^T]. In particular, in both of these 
examples, the fields of most interest are the twist fields associated to the Z2 symmetry of the 
Majorana theory ^21 These fields are the scaling limit of the spin operator in the Ising 
chain (the Pauli matrix that enjoys a nearest-neighbor interaction) and of the spin variables in 
the statistical Ising model. There are two of these twist fields: the order field and the disorder 
field. They can be seen as representing the scaling limit of the same operator, but in different 
regimes (ordered and disordered) away from the critical point. When the mass is set to zero 
(at criticality), the theory reduces to the well-known Ising conformal field theory. 

The correlation functions of twist fields exhibit non-trivial, "non-free" behaviors, and can be 
obtained as appropriate solutions to non-linear differential equations jHj. At finite temperature, 
or on the cylinder, these equations are partial differential equations in the coordinates on the 
cylinder |T5 | I16 | ITTj. and do not offer a very useful tool for numerically evaluating correlation 
functions, neither for analyzing their large-distance behavior (short-distance behaviors, on the 
other hand, can always be analyzed by conformal perturbation theory JD])- Hence it is worth 
studying in more detail the large- distance expansion of correlation functions of twist fields. 
This study constitutes a first step towards generalizations to interacting integrable models. In 
particular, we will present a method that parallels the zero-temperature form factor expansion. 



1.1 Correlation functions at finite temperature 

At finite temperature, correlation functions are obtained by taking a statistical average of 
quantum averages, with Boltzmann weights e~^^ where E is the energy of the quantum state: 



Tr [e-^^0{x, r)] 
Tr \e~ 



Here, L is the inverse temperature and H is the Hamiltonian. A consequence of the imaginary- 
time formalism ^Hl is the Kubo-Martin-Schwinger (KMS) identity l2Uj. 

{{0{x,T)...)), = {-iy{{0{x,r + L)...)), (1.2) 

where / is the statistics of O (it is 1 for fermionic operators and for bosinic operators), and 
where the dots (■ ■ ■) represent local fields at time r and at positions different from x. Then, 
finite-temperature correlation functions can be seen as vacuum expectation values in the same 
model of QFT, but this time quantized on a circle of circumference L: 

{{0{x, t) ■ ■ = e*"^/2 L(vac|OL(-r, x) ■ ■ ■ |vac)i (1.3) 

where s is the spin of O. The operator Ol{~t,x) is the corresponding operator acting on the 
Hilbert space on the circle, with space variable — r (parameterizing the circle of circumference 
L) and Euclidean time variable x (going along the cylinder). The vector Ivac)^ is the vacuum 
in the Hilbert space on the circle. There are two sectors in the quantization on the circle: 
Neveu-Schwartz (NS) and Ramond (R), where the fermion fields are anti-periodic and periodic, 
respectively, around the circle. The trace (jl.lll with insertion of operators that are local with 
respect to the fermion fields naturally corresponds to the NS sector due to the KMS identity. 
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With insertion of twist fields, however, the sector may be changed to R or to a mixed sector by 
an appropriate choice of the associated branch cuts. The phase factor in ()1.3p is in agreement 
with the fact that the operator Ol is Hermitian on the Hilbert space on the circle, if the 
corresponding operator O is Hermitian on the Hilbert space on the line. 



1.2 Previous works 

A straightforward application of the zero-temperature ideas would suggest a form factor ex- 
pansion of the vacuum expectation values on the circle ()1.3|) : 

i(vac|Oi(x,r)Oi(0,0)|vac)i = 

y , y , L(vac|OL(0,0)|ni, . . . ,nfc)i ^(ni, . . . ,nfc|OL(0,0)|vac)L 

fc=0 ni,...,nfc 

(1.4) 

where the eigenstates of the momentum operator and of the Hamiltonian on the circle are 
parametrized by discrete variables n^'s. This form is valid for any integrable models on the 
circle; it can also obviously be generalized to other multi-point correlation functions. In general 
integrable models, the energy levels £^ni,...,nfc are not known in closed form (although exact 
methods exist to obtain non-linear integral equations that define them: from thermodynamic 
Bethe ansatz techniques, from calculations a la Destri-de Vega and from the so-called BLZ 
program). Also the matrix elements of local fields l{v^c\Ol{0, 0)|?t,i, . . . , n^.) (form factors on 
the circle) still seem to be far from reach (however, see fn\ \'2'2\ 1^ )- Hence, this method 
does not seem to be applicable in general yet. Nevertheless, in the Majorana fermion theory 
with mass m, the energy levels are simply given by 



where rij G Z + | in the Neveu-Schwartz (NS) sector and rij G Z in the Ramond (R) sector (the 
sectors where the fermion field have anti-periodic and periodic conditions around the cylinder, 
respectively). Matrix elements of the primary order and disorder fields were calculated in [25] , 
in the lattice Ising model, and in a simpler way in [2Z| directly in the Majorana field theory 
using the free-fermion equations of motion and the "doubling trick" . The form factor expansion 
on the circle can then be used to obtain the leading exponentially decreasing behavior of static 
correlation functions in the quantum Ising chain, and, with slightly more work, to obtain the 
large-time dynamical correlation functions in the semi-classical regime |2H1- This reproduces 
early results of .29^ and [SU] . 

A second set of ideas was suggested by Leclair, Lesage, Sachdev and Saleur in starting 
from the representation By partly performing the trace and by using general properties 

of matrix elements of local fields, the following representation for two-point functions in the 
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Major ana theory of mass m was obtained: 



EE/ ;'n' -.^^co..., ...(".. ■ ■ ■ . ''^(^■"' ...C. .... 

(1.5) 

where F^'^ ^^(6'i, . . . , ^fc) are matrix elements of O in the Hilbert space on the hne of the form 
{\a.c\0\6j^, . . . , 9j^^, + ire, . . . ,9i-^ + in), where all 9j's correspond to ej = (+), and all 6'j's 
correspond to ej = (— ). 

This expression is better adapted to study dynamical correlation functions in real time t = it 
than the form factor expansion on the circle. Another advantage of this idea is that, following 
methods similar to those used in |22] at zero temperature, one can obtain a Fredholm determi- 
nant representation for finite-temperature two-point functions, which can be useful for further 
analysis. Such a representation was obtained in and was used in order to obtain a set of 
non-linear partial differential equations for the two-point function of twist fields. This lead to 
the known equations plus, in particular, an additional equation involving temperature deriva- 
tives. However, this additional equation does not follow from any of the other methods used in 
other works [T31 UHl Ej , and it was suggested not to hold in ^ l_7j . 

This idea also seems to allow more easily generalizations to interacting integrable models. Such 
a generalization was conjectured by Leclair and Mussardo in [33 . Saleur, however, pointed 
out that this generalization might not be correct [31] , and the numerical calculation of Castro- 
Alvaredo and Fring in the scaling Lee- Yang model tended to agree with this incorrectness |35j . 

Note that verifications of the results of [SI] were made for some fields in the complex (Dirac) 
free fermion and in the Federbush model in [33], and that agreement was found with known 
results in various limits. Leclair and Mussardo also proposed in [3S], along similar lines, a 
formula for one-point functions of local operators in interacting integrable models, which, after 
some controversy, seem to be confirmed j3ti[ I77| , although it is fair to say that more checks 
would be appropriate (on the subject of one-point functions from form factors, see also 39j). 

In the works mentioned above, very little was said concerning finite-temperature correlation 
functions with more than two operators. 



1.3 Our work 

In the following, we will develop a scheme for obtaining large- distance representations of finite- 
temperature correlation functions in the Majorana theory of the form: 



{{O,{x,T)Oli0,0))) 



L 



k=0 ei,...,efe=± ■ 



^! UU (1 + ^" 

(1.6) 
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with appropriate generically L-dependent functions ^ {9i, . . . ,9k; L). We will see that this 
formula is valid: 

• with both Oi, O2 fields that are local with respect to the fermion fields, in the NS sector; 

• with one of Oi or O2 a twist field, in a mixed R-NS or NS-R sector; 

• with both Oi, O2 twist fields, in the R sector. 

This scheme does not say anything new about one-point functions; in particular, the functions 
are evaluated up to normalization using analytic properties in rapidity space. We will only 
analyze two-point correlation functions, but the techniques and ideas introduced can be easily 
generalized to multi-point correlation functions. 

Although our expansion for two-point functions is of the same form as the expansion 

fjl.5|) . we will find that our methods for obtaining ()1.6|) is clearly more powerful than an explicit 
evaluation of the trace as in [3T], and more importantly that most of the results of (where 
this form first appeared) are in fact incorrect^. It is a simple matter to check that the positive 
verifications of the results of |3I] made in the works described in the previous subsection, being 
only for simple fields or in particular limits, can still be made on the formulas that we propose. 
We will not do that explicitly here, since we have a much stronger verification of the correctness 
of our results through their relation with the uncontroversial form factors on the circle. Our 
main results can be grouped into 4 points that exhibit the differences between our approach 
and results, and those of previous works. They are: 

1. Definition of finite-temperature form factors. The objects ^ {9i, . . . ,9k', L) in- 
volved in ()1.6|) . which we will call "finite-temperature form factors", are not given by simple 
matrix elements of O in the Hilbert space on the line as was suggested in the proposal 
PT] . They are rather given by appropriate traces of multiple commutators ([■,■]) and anti- 
commutators ({■, ■}) on the Hilbert space of the Majorana theory on the line: 

fZ...,JO„ ...,9k;L) = m^{9,), [a-(^2), {• • ■ , O{0, 0) ■ ■ -j]}))^ {9, j) (1.7) 

where a'^{9) = a'^{9) is a creation operator for a particle of rapidity 9, and a^{9) = a{9) is 
an annihilation operator for a particle of rapidity 9. Simple calculations with normal-ordered 
products of fermion fields, for instance, make it clear that generically, these objects are not 
equal to zero-temperature form factors. The difference will be explained in Point 3 below. 

The techniques we use to derive such a representation are related to the physical theory of 
particle and hole excitations above a "thermal vacuum", which was initially proposed more 
than thirty years ago and which developed into a mature theory under the name of thermo- 
field dynamics [101 (for a review, see for instance |12])- We will explain how the traces 
p.7|) can be seen as matrix elements between a thermal vacuum (on the left) and a state with 
particle and hole excitations (on the right), and how the representation comes from a 

resolution of the identity on the space of NS-sector thermal states. In fact, there are still very 

^Note that in the latter work, some singularities were neglected without good reasons (as was admitted 
there) in the derivation of H1.5|l . 
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few applications of the ideas of integrable quantum field theory (like the form factor program) 
to thermo-field dynamics (see, however, the recent study [IH] of bosonization in thermo-field 
dynamics). Our study constitutes a step in this direction. 

2. Relation between form factors on the circle and finite-temperature form fac- 
tors. We will show that an appropriate analytical continuation in rapidity space of the finite- 
temperature form factors reproduces form factors on the circle. This links the two sets of 
ideas mentioned above, and, as we will see below, gives an "analytical" way of calculating form 
factors on the circle. The precise correspondence is given by 

L(ni,...,n;|Oi(0)|ni,...,n,)i = e-^ (Y{ / ,\ ^ I I TT / ^ 1 x 

;i.8) 

where there are k positive charges and / negative charges in the indices of , and where an 
are defined by 

27rn , „ 1 , , , 

smha„ = — - {neZ + -). (1.9) 

TTXlj A 

This formula is valid for any excited states in the NS sector. When O is a twist field, formula 
()1.8|1 can then be applied only if one of the bra or the ket is the vacuum, and if the branch 
associated to this twist field is chosen so that this vacuum is in the R sector and the excited 
state is in the NS sector. Using this formula, we will show that the form factor expansion on 
the circle and the finite-temperature form factor expansion give the same result. 

3. Evaluation of finite-temperature form factors. We will evaluate explicitly the finite- 
temperature form factors f^^,,,^^^^!^ ■ ■ ■ ,(^k]L) for some fields. We find: 

1. Mixing. For any local non- interacting field O (a field that is local with respect to the 
fundamental fermion field), they are equal to matrix elements in the Hilbert space on 
the line of a sum of a finite number of local fields, including O as well as operators of 
lower scaling dimensions and of equal or lower conformal dimensions. We will present 
a procedure for calculating the associated mixing matrix. This mixing can ultimately 
be seen as coming from the difference between the normal-ordering on the line and the 
normal-ordering on the circle. Note that in the proposal f3I], only the energy field, among 
non-interacting fields, was studied, and the mixing was treated by adding to this field the 
identity operator in such a way that it have zero thermal expectation value. For other 
fields, such a simple prescription would not be enough. 

2. Leg factors. For the order and disorder fields, they are given by their matrix elements in 
the Hilbert space on the line times functions depending on the individual rapidities ( "leg 
factors"). These factors were absent in |31j . 

Note that, consistently with ()1.8|) . both the mixing phenomenon and the appearance of leg 
factors are observed in form factors on the circle. In particular, we were able to reproduce from 
our results the known form factors on the circle. 
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We evaluate the finite-temperature form factors of order and disorder fields by deriving and 
solving a Riemann-Hilbert problem that completely fixes the set of finite-temperature form 
factors of all local twist fields (order and disorder fields and their descendants) in the Majorana 
field theory. The Riemann-Hilbert problem depends explicitly on the temperature. For descen- 
dant of order and disorder fields, we present a conjecture for the structure in which the mixing 
phenomenon can be described (without, though, calculating the associated mixing matrix). 

It is important to note that since the effect of temperature on the finite-temperature form factor 
of order and disorder fields is only in the leg factors, a Fredholm determinant representation 
can straightforwardly be obtained from the finite-temperature form factor expansion, as it was 
done in |31j. One needs only replace the filling fraction of [31j by the filling fraction times our 
leg factors. However, because of these factors, following the derivation of jSI] it is clear that 
the equation involving temperature derivatives is incorrect, in accordance with the suggestion 



4. Resolution of the divergencies. Because of the poles in finite-temperature form factors 
of twist fields, the representation ()1.6|) . as it stands, suffers from divergencies at colliding 
rapidities, when they are associated to opposite values of e. Using the link, provided by ()1.8|1 . 
between this representation and the form factor expansion on the circle, we will describe the 
correct prescription that makes the integrals in ()1.6|) finite. This point was left slightly unclear 
in the proposal [31j, were a similar problem arose. Such divergencies at colliding rapidities, 
however, were studied more carefully, in a related context, in the recent paper j2H] using ideas 
similar to those presented here. 

Some of the ideas that we will present have an easy generalization to interacting integrable 
models; this raises hopes that finite-temperature form factors can also be evaluated in such 
models. Work in this direction is in progress. 

The paper is organized as follows. In Section 2, we recall the quantization of the Majorana 
theory on the line and on the circle. In Section 3, we introduce the space of "finite-temperature 
state" and explain how finite-temperature form factors give large-distance expansions of corre- 
lation functions. In Section 4, we give an argument from basic principles justifying the relation 
between finite-temperature form factors and form factors on the circle. In Section 5, we study 
finite-temperature form factors of non- interacting fields. In Section 6, we derive the Riemann- 
Hilbert problem for finite-temperature form factors of twist fields. In Section 7, we calculate 
these finite-temperature form factors. Finally, we conclude with a series of open problems. 

2 Quantization of the Majorana theory on the hne and 
on the circle 

In the free Majorana theory with mass m quantized on the line, fermion operators evolved in 
Euclidean time r are given by: 



of nzi. 
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where the mode operators a{6) and their Hermitian conjugate a'^{6) satisfy the canonical anti- 
commutation relations 

{a\e),a{e')} = 5{e-e') (2.2) 

(other anti-commutators vanishing) and where 

pe = msinh^ , 
Eq = m cosh 6 . 

The fermion operators satisfy the equations of motion 

— 1 TTl - 

dip{x, r) = - (9x. + idr)i' = —ip 

— 1 — HI 

di,{x,T) ^ -{d,-rd^)^ = -i, (2.3) 
and obey the equal-time anti-commutation relations 

{z^(x), ^(x')} = 5{x - x') , {V^(x), V^(x )} = 5{x - x) , (2.4) 
which is simple to derive from the representation 

8{x) = dpe'P'' (2.5) 

of the delta-function. The Hilbert space Ti. is simply the Fock space over the algebra ()2.2|1 with 
vacuum vector |vac) defined by a(6')|vac) = 0. Vectors in H will be denoted by 

\9^,...,9k) = a\e,)■■■a\9k)\Y^c) . (2.6) 

A basis is formed by taking, for instance, 6i > ■ ■ ■ > 6k- The Hamiltonian is given by 

/oo 
dOm cosh ea^{e)a{e) (2.7) 
-oo 

and has the property of being bounded from below on Ti. and of generating time translations: 

[H, ij{x, t)] = —i:{x, t) , [H, i,{x, r)] = — ^(x, r) . (2.8) 

For future reference, note that the leading terms of the OPE's ipi^: t)'V^(0, 0) and ip{x, r)-?/'(0, 0) 
are given by 

^(x, r)V'(0, 0) ~ _ . ' . , , ^(x, r)V'(0, 0) ~ / . . . (2.9) 

This normalization implies that the Majorana fermions are real, and is different from the usual 
normalization used in conformal field theory. 
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On the other hand, in the same theory quantized on the circle of circumference L, with anti- 
periodic conditions on the fermion fields, the fermion operators evolved in Euclidean time r 
are: 



Mx,t) = --L= J2 (an e-^"-^"- - at e-^"-+^"-) (2.10) 

where the discrete mode operators a„ and their Hermitian conjugate a'l satisfy the canonical 
ant i- commutation relations 

{al, ttn'} = 5n,n' (2.11) 

(other anti-commutators vanishing) and where 

27m ,1, 
Pn = msinha„ = — — (nGZ+-), (2.12) 

Ij Z 

En = mcosha„ . 

The fermion operators satisfy the equations of motion ()2.3p as well as the equal-time anti- 
commutation relations ()2.4j) (with the replacement ip ^ ipi and ip ^ ^Pl)'-, the latter is simple 
to derive from the representation 

5{x) = \Y^ e^''"" (2.13) 

of the delta-function, valid on the space of antiperiodic functions on an interval of length L. 
The Hilbert space is simply the Fock space over the algebra ()2.11|) with vacuum vector 
|vac)/, defined by a^lvac)^ = 0. Vectors in Tii will be denoted by 

|ni,...,nfc)L = 4^---4jvac)L . (2.14) 

A basis is formed by taking, for instance, ni > ■ ■ ■ > n^. The Hamiltonian is given by 

Hl = mcosha„ata„ (2-15) 

and has the property of being bounded from below on Tii and of generating time translations: 

d - d - 

[Hl,iI)l{x,t)] = —iPl{x,t) , [HL,ipLix,T)] = — V^i(x,r) . (2.16) 

The leading terms of the OPE ip^^x, T)ip ^{0,0) and ijjL{x,T)ijjL{0,0) are the same as in ()2.9p . 
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3 The space of "finite-temperature states" and large- 
distance expansion of correlation functions 

We are interested in calculating correlation functions in the theory with Hamiltonian H at 
finite temperature 1/L. They are given by ratios of traces on Ti as in (jl.ll) . In particular, 
the relation (jl.Hj) holds, where for non-interacting fields (normal-ordered products of fermion 
fields and of their space derivatives), the operator Ol can be defined by first writing the 
regularization of O in terms of a product of fermion operators ip, ip at different points (without 
normal ordering), then by replacing all fermion operators by their corresponding forms on the 
circle: ^ i^L, i' ^ V'l- 

The space of "finite-temperature states" that we will use is in fact well-known in the literature: 
it is the Liouville space jS], or the space of operators of the theory, End(H). As recalled in the 
introduction, the study of finite-temperature quantum field theory using the Liouville space 
goes under the name of thermo-field dynamics. It is not our intention to go into any detail of 
this field of study; we will recall only basic principles, in a slightly different, but equivalent, 
formulation, then develop certain parts for our purposes. 

We will denote by C the (completion of the) infinite-dimensional subspace of End('H) spanned 
by products of any finite number of mode operators a{6), a^{6) at different rapidities. This 
subspace contains local non-interacting fields of the Majorana theory. The vacuum in C is 
defined by the identity on H: 

|vac)^ = In , (3.1) 
and a complete basis of states is given by 

k 

K ■ ■ ■,^k)t...,e, = n + e"'''''''^) ■ ■■a^'i^k) , e,>e2>--->ek (3.2) 

where ej are signs (±), which will be called "charge" of the "particle" of rapidity 6j. In fact, 
the sign ej is associated to the creation from the thermal bath, if it is positive, or to the 
absorption by the thermal bath, if it is negative, of a particle of rapidity 6j. Note that in (|3.2p . 

we introduced the factor 11^=1 + e~^^'"^'^i j ; this just specifies the normalization of the states, 
but will play an important role in the following. The inner product on C is defined by 

{u\v) = {{UW))l , for \u) =U, \v) = V. (3.3) 

In fact, as usual, it will be convenient to define states with other orderings of the rapidities by 

1^ nxC _ / (-l)^l^^'{l)'---'^P(fc))fp(i),...,.p(,) (^iT^^jVz^j) /o4^ 



i3t ^j\e, = e^ 



3/ 



Here P is a permutation of the k first integers such that Op{i) > ■ ■ ■ > Op^j.-^, and (—1)'^ is the 
sign of the permutation. Operators on Ti can also be mapped to operators on £ by implementing 
their action on C through their left-action or their right-action on endomorphisms of H.. It will 
be convenient for our purposes to concentrate solely on their left-action: 

O e End(H) e End(£) = left-action of O on End(H) . (3.5) 
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Using the cyclic property of traces, it is easy to calculate the following trace: 



{{a{e)a\9')))L = e^^«((at(0')«W))i 



{-{{a{e)a\e'))), + s{e-e')) 

5(9 - 9') 



'LEo 



(3.6) 



Similarly, 



and we also have 



{{a'{0H9')))L = (3.7) 



{{a{9)a{9'))), = {{a\9)a\9'))), = 0. (3.8) 

From these results and from using Wick's theorem (which applies for traces of products of free 
modes), it is easily seen that 

k 

,,_i{9,, ...,9,\9[,..., 9i)f,_,^ = 5{9, -9[)... S{9, - 9',) 5,,,. ■ ■ ■ <5,,,. J] (l + ^"^^'''O 

i=i 

(3.9) 

if > ■ ■ ■ > ^^fc and 9[ > ■ ■ ■ > 9'i^, and that other inner products with such an ordering of 
rapidities are zero. 

From these definitions, we see that finite-temperature expectation values of operators on Ti are 
vacuum expectation values on C: 

{{0{x,t)---))l= ^(vac|C^(x,r)---|vac)^ . (3.10) 

Hence, using the resolution of the identity on C: 

^c = ^Yl [ 7777r~7 ^^-eLE. \ , (3.11) 



fc=0 ei,...,efc 



two-point functions can be expanded as in ()1.6p . where we define finite-temperature form factors 
as 

fZ...,Jdu---,dk;L)= ^(vac|O^(0, 0)1^1,..., e,)f^_,^ . (3.12) 

Of course, multi-point correlation functions can be expanded in similar ways. Note that in 
obtaining (jl.6|) . we used the symmetries under translations (along the cylinder and around the 
cylinder) in order to bring out the x and r dependence of 0{x, r) in exponential factors. 

It is easy to check, using the cyclic property of traces, that the finite-temperature form factors 
can be written in the form (jl.7|) . This will be useful in understanding why these matrix elements 
are related to form factors on the circle. 

It is a simple matter to observe that in the zero-temperature limit L —>■ oo, the finite- 
temperature form factors reproduce the zero-temperature form factors: 

lim f^^ _{9i, . . . , 9^+, 9k++i, ...,9k;L) = {9k, . . . , 9k+^i\O{0, 0)\9i, . . . , 6^^^) 

^-^"^ (3.13) 

{9,^9,\/te{l,...,k+},je{k^ + l,...,N}) 
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where there are fc+ positive charges (+)'s in the indices of f^. For other orderings of the charges 
(and of the associated rapidities), we get extra minus signs. The finite-temperature form factor 
expansion ()1.6p . in the hmit L —>■ oo, then reproduces the usual form factor expansion, since 
only the form factors with ej = + for all j's remain. 

It is also a simple matter to observe that more general matrix elements in the space C are 
simply related to the finite-temperature form factors ()3.12|) : 

fZ...,eS^l^ • • • ) (^k+,dk++l, ■■■,Ok]L) = ^efc,.. -6fc_^+^i(6'A:, . . . , 6'fc++i | O'^ (0, 0)|6'i, . . . , )fj 

(^,^^,V2G{l,...,A;+},jG{A;+ + l,...,/c}) . 

(3.14) 

Let us make two comments. First, note that expression ()l.(i|l . as it stands, is not a practical 
large-distance expansion. It is more convenient, as in the zero-temperature case, to analyti- 
cally continue to complex values of ^/s, in such a way that the argument of the exponential, 
i X]j=i ^jPdjX, acquires a real part that becomes large and negative at large \6j \ for any j. This 
is optimal, if x > 0, at 9j ^ 6j + iej'K/2. However, in contrast to the zero-temperature case, 
the theory is not invariant under rotation, hence this shift of rapidities does not keep the form 
factors invariant. We will see how this shift is related to a form factor expansion on the circle. 

Second, the integrals in the expansion ()1.6p may suffer from singularities at colliding rapidities 
in their integrands. As we will see, such singularities arise in finite-temperature form factors 
of interacting fields when rapidity variables are associated to opposite charges. We need a 
prescription to define properly the finite-temperature form factor expansion in such a case. The 
correct prescription is the following: Consider the difference x between the positions of the first 
and of the second operator involved to be positive; then, we simply shift slightly the integration 
lines towards the positive or negative imaginary direction in rapidity planes, in such a way that 
the integration stays convergent at r = and that the poles are avoided. More precisely, the 
integrals on 6'i, . . . , 9k in fll.6|) should be on contours parallel to the real line with Im(6'j) = ejO+. 
We verify that this prescription is correct, by showing independently, in Appendices 1X1 and iBl 
the formula ()1.8|) . and the fact that our prescription makes the finite-temperature form factor 
expansion (jl.6|) and the form factor expansion on the circle (jl.4p equivalent. 

A note on Ramond and Neveu-Schwartz sectors. 

Because of the relation (|1.2j) . fermion operators satisfy anti-periodic conditions in Euclidean 
time under the trace. The trace (jl.ip with insertions of fermion fields then corresponds to 
correlation functions on the cylinder with anti-periodic conditions for the fermion fields around 
the cylinder, that is, in the NS sector in the language of the quantization on the circle. 

When twist fields are inserted inside the trace, the situation is slightly different. Twist fields 
are end-points of branch cuts (or line defects), and on the cylinder, a semi-infinite branch cut 
can be in two directions. The definition of operators corresponding to twist fields is recalled 
in Appendix O we denote two types of operators corresponding to any twist field O by (9+ 
and with branch cut towards the positive and negative x direction along the cylinder, 
respectively. We call these operators "right-twist" and "left-twist" operators, respectively. 

Finite-temperature correlation functions (|1.1|) with insertion of only one twist fields (and any 
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number of fermion fields) correspond then to correlation functions on the cylinder in a mixed 
sector, the precise mixed sector depending on which operator (right-twist or left-twist) is taken. 
The fermion fields have periodic conditions around the cylinder (R sector) where the branch cut 
stands, and anti-periodic conditions (NS sector) elsewhere. This property is essential in deriving 
the Riemann-Hilbert problem for finite-temperature form factors in Section IHl Correspondingly, 
our formula (jl.fij) can be used with any one of Oi and O2 a twist field, and in both cases, the 
branch cut can be in either direction. 

Finite-temperature correlation functions with insertion of two twist fields can correspond to 
correlation functions on the cylinder in a pure NS sector or in a pure R sector, again depending 
on which operators are taken under the trace. Our formula ()1.6|) . however, can only be used, in 
this case, to give an expansion of correlation functions in the R sector. This can be interpreted 
by the fact that in each finite-temperature matrix element of a twist field, there is no other 
twist field to provide another end-point of the branch cut, so that it must go to infinity along 
the cylinder, hence it changes the natural NS sector of the trace to an R sector at infinity. 
This can be seen technically by analyzing our Riemann-Hilbert problem given in Section IHl and 
by doing a contour-deformation analysis of the expansion ()1.6p . similarly to what is done in 
Appendix |B| Correspondingly, then, the expansion ()1.6|) is valid, if both Oi and O2 are twist 
fields, only when they have branch cuts in opposite directions. 

In the same spirit, note that formula ()1.8|) for calculating form factors on the circle from finite- 
temperature form factors, explained in the next section, can only be used with the vacuum (in 
the quantization on the circle) in the R sector, and excited states in the NS sector. 

However, it can be useful to note that one could "twist" the construction of this section by 
defining the inner product as traces with periodic conditions on the fermion fields: 

where the operator 0"+ is the operator associated to a primary twist field of non-zero finite- 
temperature average, with a branch cut on the right. The resulting finite-temperature corre- 
lation functions of fermion fields will be in the R sector. With such a definition, the formula 

fll.7|) is still valid, but the normalization factors 11^=1 (^1 + 6'^^^^"^^ introduced above should 

be replaced by 11^=1 (^^ ~ e"^^^^*-' and a similar replacement should be made in the finite- 
temperature form factor expansion Following a similar reasoning as above, the resulting 
expansion for the two-point function of twist fields will be valid, in this construction, only in 
the NS sector. The finite-temperature form factors obtained in this construction can then be 
used to calculate form factors on the circle with excited states in the R sector and vacuum 
in the NS sector, with a formula similar to ()1.8|) . Note that starting from the form factor 
expansion on the circle, a representation similar to (more precisely, after some approxi- 
mations, to (jl.5|) ) was obtained in |2.8] for two-point functions of order fields, with the factors 

Y[^=i + e'""^^^"^^ replaced by 11^=1 (^1 ~ e^'^^'^^'^j^ . The construction just outlined in this 
paragraph explains this phenomenon. In the rest of the paper, though, we will not consider 
this construction. 
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4 From finite-temperature form factors to form factors 
on the circle 



4.1 General idea 

Observe that the two equivalent expressions (jl.lll and ()1.3|) for correlation functions at finite 
temperature, or on a cylinder, result from two different quantization schemes for the same 
theory: in the first, the equal-time slices are lines along the cylinder, whereas in the second, 
they are circles around it. 

Observe also that the finite-temperature form factor expansion of two-point functions (jl.fij) 
could be made into a form factor expansion on the circle ()1.4j) . if finite-temperature form 
factors have appropriate properties, by deforming the contours in rapidity space and taking the 
residues of the poles of the measure. 

There are more precise observations that allow us to infer that form factors on the circle can 
be obtained from finite-temperature form factors. Let us develop two of them below. 

First, consider a massive free theory on the Poincare disk (the Dirac theory on the Poincare 
disk was studied in |45j, and we present here ideas exposed there), and consider a quantization 
scheme where equal-time slices are the orbits of one of its non-compact isometry (that is, the 
momentum operator generates isometry transformations X and its spectrum is continuous). 
In this quantization scheme, consider the set of states |P) that diagonalize the momentum 
operator, with real continuous eigenvalues P. In fact, the states can be described by a set of 
momenta {p} and a set of discrete variables (particle momenta and particle types). 

Now consider matrix elements of local spinless operators in this set of states, of the form 
(vac|(9|{p}), and see them as analytically continued functions of momenta {p}. Choose a 
normalization of the states in such a way that these matrix elements are entire functions of the 
momenta {p}'s. With this choice, there is a measure involved in the resolution of the identity: 
1 = l\^V\ p({p})|{p})({p}| where k is the number of particles. The singularity structure of 
p({p}) then describes the energy spectrum of the theory in the quantization scheme where the 
roles of the momentum and the Hamiltonian are exchanged: the isometry X is generated by the 
Hamiltonian. In particular, poles at purely imaginary values p = -^lE correspond to single- 
particle discrete eigenvalues E of the Hamiltonian in that scheme. Also, the matrix elements 
of O in the "exchanged" quantization scheme, of the form '(vac|(9|{i?})' where Ij-E})' is an 
energy eigenvalue, are proportional to A/Res(p({p})) (vac|(9|{p}) at the values {p = ±iE}. 

The second observation is a similar one, but this time for interacting relativistic models on flat 
space. Consider for simplicity such a model with only one particle type (one mass). Consider a 
set of momentum eigenstates, parametrized by rapidity variables {6}. Choose a normalization 
of the states such that the form factors of any local field are entire functions of the total rapidity 
when rapidity differences are fixed, and such that the measure involved in the resolution of the 
identity 1 = "^i. J[d0]p{{9})\{9}) {{9}\, where k is the number of particles, factorizes into a 
product of measures for the individual rapidities of each particle. This is the usual choice: the 
dependence of form factors of local fields on the total rapidity is the entire function e^^^^\ 
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where s is the spin of the field; and the measure is p{{0}) = 1. To clarify the idea, the resolution 
of the identity can be written as integration over momentum variables with an appropriate 

dp 



relativistically invariant measure: 1 = J 



The measure has branch 



cuts from p = ±im towards ±ioo: these describe, again, the energy spectrum of the theory 
where the roles of the momentum and of the Hamiltonian are exchanged. Of course here, by 
relativistic invariance, this exchange gives the same Hilbert space. It can be implemented by 
the shift 6 6 + iTT/2 of rapidity variables, and the form factors are invariant, up to a phase 
due to the spin, under such a shift of rapidities. 

Similarly, it is not hard to relate analytical continuation of special traces in angular quantiza- 
tion to the spectrum and matrix elements of local operators in radial quantization. Adapting 
these observations to the case of the cylinder, the main idea for obtaining form factors on the 
circle from finite-temperature form factors in the free Majorana theory can be expressed in the 
following (schematic) steps: 



• Consider the Liouville space £, where the trace 

Tr (e-^^O(s,0)C(0,0)) 
Tr (e-^^) 

is a vacuum expectation value and where eigenvalues of the momentum operator are 
described by continuous variables, the rapidities 9j; 

• Find a measure p{{9}) in 

k=0 {e} ^ 

such that matrix elements of local fields ^ which are also local with respect to the fun- 
damental field Tp, ■^(vac|\l/(0, 0)|{6'})|;|, are entire functions of the rapidities {6*}; 

• Calculate form factors on the circle by analytical continuation in the rapidity variables 
to the positions a„ + 27r/2 of the poles of the measure p: 

L(vac|^L(0,0)|{n})LOC VResp ^(vac|^(0, 0)|{a„ ± m/2})f . . 



In the last step, a„ are defined in ()1.9|) . In the second step, the requirement that the fields 
\l/ be local with respect to the fundamental field ip comes from the fact that rotation is not a 
symmetry on the cylinder. Finite-temperature form factors of twist fields (which are essentially 
the end-point of linear defects) can be expected to have a more complicated pole structure. 

More than relating finite-temperature form factors to form factors on the circle, the steps above 
tell us that an appropriate choice of the measure is related to a simple analytical structure of 
finite-temperature form factors. 

It is hopefully possible to generalize some of these requirements to interacting integrable models, 
but this will not be pursued further here. 
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4.2 Implementation of the general idea 



For the massive Majorana theory, the first step was described in the previous section. 

We now show that the second step is reahzed by the definition ()3.12j) . due to the representation 
()1.7|1 . Consider a complete basis j} of local operators that are also local with respect to the 
fundamental field ip] these are non-interacting operators, for instance normal-ordered products 
of fermion operators and of their space derivatives, with coefficients that may contain integer 
powers of the mass m (the mass m has scaling dimension 1 and spin 0). We consider the basis 
j} to be for the linear space of operators on the field of complex numbers without dependence 
on m. We can take the operators j to have well-defined spins Sj and scaling dimensions dj. 
With the fermion operators iIj and ip, they satisfy simple equal-time commutation (or anti- 
commutation) relations: 

[^(x),vl/,(x')] = 

j; dj=0,...,di-| 

[V'(x),^i(x')] = 4 m ,{x') 6^''^-''^--^\x - x') . (4.1) 

j; dj=0,...,di-| 

Here the symbol [■, ■] is the commutator if Sj is an integer or the anti-commutator if Sj is a half- 
integer (the spin Sj determines the statistics of the operator j{x')). The sum is of course finite, 
because there are no fields of negative dimension. The coefficients c^, do not depend on m. 
The fact that ^ j have well-defined spin and dimension imposes constraints on the coefficients 
c^, c^, but we will not use these constraints here. 

The operators a{0) and a^{6) can be decomposed in terms of the local operators tpi^x) and ■?/'(x): 
a\e) = -.f^ r rfxe^P"" (e^/2^(x) - tie-'l^i){x)) . (4.2) 



2 V „ _ 

The one-particle finite-temperature form factor of a field Oi{x) can then be written, from ()1.7|) . 
ff^{e-L) = dxe"^^- {e'l\m^).'^^m))L~e^e-'/\m^).'^^m))L) 

1 trf '■^ 



2V vr ^ 

di 

1 m 



2S E i-^^Por^-''-^' - "e-^/^c^) {{^m))L (4.3) 

Note that these finite-temperature form factors are entire functions of the variable 6. 
Similarly, the multi-particle form-factors are given by 

/*:':'.(9i,...,9i;i)= (4.4) 

^ ^ ^ ji,-,jk 1=1 
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Again, they are entire functions of the variables 6i. Hence, the measure appearing in the 
resolution of the identity ()3.11|) fulfills the condition of the second step. 

The third step then tells us how to evaluate form factors on the circle. The precise relation can 
be obtained, for instance, by calculating explicit examples for simple non-interacting fields. It 
is given by ()1.8|) . A more direct derivation of this relation is presented in Appendix EI A proof 
that the finite-temperature form factor expansion of two-point functions ()1.6|) is equivalent to 
the expansion in form factors on the circle ()1.4|) is presented in Appendix ^ using properties 
of finite-temperature form factors that are derived in the rest of this paper. 



5 Finite-temperature form factors of non-interacting fields 
5.1 General properties and particular cases 

The expression fl4.4|l gives a general formula for finite-temperature form factors of local non- 
interacting operators. In particular, it is an easy matter to verify the property 

/*:...,.,,...,.. (^1' ...,% + ^vr, 4; L) = tfZ..,-.„...,e,iOu . . . , . . . , 4; L) . (5.1) 

This leads to the quasi-periodicity property 

• • • ' + 2*^' ■■■.Ok;L) = -/*:...,,^,...,,,(^i, ...,e,,...,e,;L) . (5.2) 

Note that the same formula ()4.4|1 could be used, with the replacement {{'^j))L ^ (^i); for 
calculating form factors at zero temperature; indeed, the parameter L is only involved in 
the expectation values ((\E'j))i. This suggests that finite-temperature form factors of non- 
interacting operators are "not so far" from their zero-temperature limit. For instance, the 
finite-temperature form factors of the fundamental fermion operators are 

ftiO;L) = \^e'l^ , fi{e-L) = T^\^j^e-''- , (5.3) 

which are temperature independent. On the other hand, for the energy field (from the viewpoint 
of the Ising field theory) e = i : ipip :, the two-particle finite-temperature form factors are 

n^^{9^,92;L) = -r_,_{e,,e,-L) = £sinh (5.4) 

and 

/;^_(0i, e^- L) = -r_^4e„ e,; L) = cosh (^^) (5.5) 



which agree with the zero-temperature form factors, but the thermal expectation value is 

m dO 
TT 7n 1 + e 



n-.L) = {{e)), = - I (5.6) 
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which is non-zero and L- dependent. Those in fact are the zero-temperature form factors of the 
field e + {{e))i,l. One could re-define the field e by subtracting this thermal expectation value, 
but for a more general non-interacting field, such a subtraction will not bring all its finite- 
temperature form factors equal to its zero-temperature form factors. This aspect was missing 
in ISIj. In general, fields that are normal-ordered products of k fermion fields have fc-particle 
finite-temperature form factors independent of temperature, but generically have non-zero, 
temperature-dependent j-particle form factors for j < k {j and k being of the same parity). In 
order to describe this, we introduce the concept of mixing, which is a simple generalization of 
the description e ^ e + {{e)) L'^ of finite-temperature form factors of the energy field in terms 
of its zero-temperature form factors. 

5.2 Mixing 

Consider a complete basis {"^a} of local non-interacting operators on the field of polynomials 
in m - we will use indices a,b to label all elements of this set, instead of the indices i,j 
used in the previous section to label elements of a different basis. The elements of this set of 
operators can be taken to be in one to one correspondence with a basis of non-zero operators 
in the free massless Majorana conformal field theory. A convenient basis {^E'a} is obtained by 
considering normal-ordered products of fermion operators ^p and ip and of their holomorphic or 
anti-holomorphic derivatives, d'^ip, d^ip, with coefficients that are independent of the mass m. 
We will consider their holomorphic and anti-holomorphic dimensions A^, A^, related as usual 
to their scaling dimensions rf^ = A^ + A^ and their spins = A^ — A^. 

We will show that: 



■■■.Oh. ■■■.Ok,L) = (4, . . . , ^'^~'"M,^(mL)vl/,(0)|^i, . . . , J 



where there are /c+ positive charges (+) in the indices of Here the mixing matrix M^(mL) 
mixes operators '^h of lower dimension than that of \l/a, db < da, of equal or lower associated 
conformal dimensions, A^ < A^, A^ < A^, and of the same statistics. The operator is a 
descendant under the fermion operator algebra of all operators ^^j,; in other words, the mixing 
occurs between and some of its ascendants of the same statistics. The sum over b is finite. 

One can see the space £ as a Fock space over the algebra of modes A^{6) and of their Hermitian 
conjugate on C, Al{9): 



{Alie),AAe')} = ii + e-^'^^'>)6ie-e')6,,, , {Ali9),Al,ie')} = {Aie),AAe')} = o . (5.8) 



b 



ie,j^e,Wze{l,...,k+},je{k+ + l 



.,k}) 

(5.7) 



The vacuum is defined by A^{9)\va.c)'-' = 0, and 




(5.9) 




(5.10) 
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One can immediately verify, for instance, that this representation on L gives {[a'^(6')]^, [a^'(^^')]^} = 
bid — 6') (Jg^e', as it should. 

Now consider an explicit set of local non- interacting operators {^a} given by the coefficients 
in the expansion around xi = 0, . . . , = of the generating operator 

^p(a;i, . . . = : ieiXi) ■ ■ ■ '4)p{Xp,iepXp) : (5.11) 

where each element of the set of fields {tpj, 1 < j < p} is one of or tp, and where ej = —1 if 
ipj = ip and ej = 1 if tpj = ip. According to the corresponding operator acting on C by 

left action is denoted by \E'p ^(xi, . . . , Xp+q). The normal-ordering operation : ■ : on 7i brings the 
operators a{9) to the right of the operators a^{9). This normal ordering operation can also be 
defined on operators acting on C by bringing operators Aj^{6) and AL(^) to the right of A^_^{6) 
and A_{6), so that we can write 

-^'j;{xi,...,Xp) = :'4)^{xi,ieiXi)---'ip^{Xp,iepXp): . (5.12) 

We can also define a more natural normal-ordering operation on C which brings all operators 
A+{e) and A_{9) to the right of ^4^(6^) and ^4^(6'); let's denote it by : ■ We introduce the 
generating operator 

^p{xi, ...,Xp) = : tjji{xi,ieiXi) ■ ■ ■tpp{Xp,iepXp) :£ = : ^p(a;i, ...,Xp):c ■ (5.13) 
By Wick's theorem, we have 

^ p ("^l' • • • ! '^p) ~l~ ^ ^ ( 1) CfYi,n{j^m -^n) 2 ("^l ! • • • i • • • ; i-ni • • • i "^p) 

(5.14) 

^ p ("^l' • • • ! "^p) ~l~ ^ ^ ( 1) CfYi,n{j^m "^n) ^p_2 ("^l ! • • • ? -^rrt; • • • ; '^n; • • • ; "^p) • 

Here Cm. 

,n{xm ^n) IS the vacuum expectation value on T~C of the product of the fermion field at 
Xm and thejermion field at Xn] that is, Cm,n{^m-Xn) = {yac\il)m{xm,iemXm)4^n{xn,ienXn)\ysic) . 
Similarly, Cm,n{xm — Xn) is the vacuum expectation value on C of the product of the field at 
Xm and the field at a;„; that is, Cm,n{xm - Xn) = ^{vac\ilj^{xm,iemXm)'^n{xn,ie„Xn)\vac)^ = 
{{iprnixm, i^mXm^'ipnixn, i€nXn))) L- The dots (. . .) indicate terms with more and more contrac- 
tions, until all fields are contracted, in the usual way. 

By solving iteratively ()5.14j) . we can write 

^ p ("^l' • • • 5 Xp^ ^ p ('^1' ■ ■ ■ 5 "^p) ~^ ^ ^ -^'m,n{Xm "''n)^p— 2('^l5 • • • ) Xm^ • • • 5 -^n; • • • ; -^p) 

m,n, m<n 

(5.15) 

where M.ra,n{xm - Xn) = (-1)™ " ^{Cm,n{xm " Xn) - Cm,n{xm " Xn))- The dots (. . .) represent 
terms containing operators with decreasing total number of fields. 

Consider the matrix element 

^(vac|^J(xi,...,Xp)|^i,...,^p)f^_,^ . (5.16) 
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Since there are as many rapidities in the state as there are fermion operators in the operator 
of which we take the matrix element, the normal ordering does not affect the result (there can 
be no internal contractions). Hence we have 

^(vac|^J(xi, . . . , . . . , ^p)f,,...,,^ = ^(vac|*^(xi, . . . , . . . , ^p)f,,...,,^ • (5.f 7) 

Moreover, by dimensional analysis, in the expression ()4.4j) for the right-hand side of the equation 
above, only vacuum expectation values of operators of dimension will remain: the identity 
operator. Since the finite-temperature expectation value of the identity operator is independent 
of the temperature, and since the only L-dependence in fl4.4|) occurs in expectation values, 
we can specialize the right-hand side of the previous equation to zero temperature (L — > oo) 
without change in the equation. It is convenient at this point to take, without loss of generality, 
the p+ first charges to be positive, and the rest to be negative, and to assume Oi ^ 6j 'i i & 
{!,... j E {p+ + I, . . . ,p}. Then we find: 

^(vac|^'^(a;i,...,Xp)|^i, 

Now, all matrix elements of . . . ,Xp) between the vacuum and states containing more 

or less than p rapidities are zero because of the finite-temperature normal ordering. Similarly, 
all zero-temperature matrix elements of \l/p(a;i, . . . ,Xp) between states whose total number of 
particle is more or less than p are zero. Hence we can write in general: 

^{ysic\'^p{xi, . . .,Xp)\9i, 
for any k and any k^. Using these equalities and the relation ()5.15|) . we have 

; ; • • • ; ) • • • ) -^p) | ^1 ) • • • 5 ) ~l~ 

(5.20) 

This is exactly of the form ()5.7p . once expanded in powers of x^'s. From this calculation, we 
see that the elements of the mixing matrix can be seen as coming from the difference between 
the normal-ordering : ■ :£, which for local fields on Ti. is essentially a normal-ordering on the 
circle, and the usual normal-ordering on the line : ■ :. 

It is possible to describe the mixing matrix in a more precise fashion. Recall that the space C 
contains local non-interacting fields. The mixing matrix can be seen as an operator acting on 
£, such that a state in C corresponding to the field ^E'a(O) is mapped to the state corresponding 
to ^^L^''-'^^M^{mL)<ifbiO). 

Consider the map ^E'(O) ^— ^^(0) that maps any normal-ordered non-interacting operator \E'(0) G 
End(7Y) to the operator ^'(0) G End(H) such that ^'^(0) = : ^'^(0) :£ (the finite-temperature 



. . . , Op^, 9p^-\-i, ■ ■ ■ 

{Op, . . . , ^p++i|^p(a;i, . . . , a;p)|6'i, . . . , 9p ) 



(5.18) 



• • • , Ok+-,Ok++i-, ■ ■ ■ ..,+,-,...,- — 

{Ok, • • • , Ok.+il'^pixi, . . . , Xp)\6i, . . . , 9k ) 



(5.19) 



n^p{xi,...,Xp) ,r) a a a . T\ _ 

_,,,,_(,c7i, . . . , yfe+, ■ ■ ■ i^k, -t^j — 

{Ok, ■ ■ ■ , ^fc++i|^p(a;i, . . . , Xp)\Oi, . . . , Ok^) + 

Mm,niXm " Xn){Ok, 61^^+1 1 ^p_2 (^1 , . . . 
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form factors of ^t'(O) are equal to the zero-temperature form factors of ^1/(0), as described above). 
We saw above how to construct ^{0): it indeed exists, and it is unique. Now consider the state 
= vl/^(o)|vac)^ as well as the state |^)^ = ^^(0)|vac)^ in C. The mixing matrix M^(mL) 
is defined by 

(5.21^ 



Define the operator U G End(£) generating this mixing matrix: 
From the arguments above, it can be shown that 



U = exp 



d9- 



1 + e^^e 



(5.22) 



(5.23) 



6 Riemann-Hilbert problem for finite-temperature form 
factors 

6.1 Non- interacting fields 

In Section we described finite-temperature form factors of non-interacting fields. They are 
all given by entire functions of the rapidity variables ^j's, as are zero-temperature form factors 
of non-interacting fields; the main non-trivial phenomenon is that of mixing between a field 
and its ascendants. 

In fact one can verify that the linear space of functions /* of the variables 9i, i = 1, . . . , k with 
the properties: 

1. acquires a sign under exchange of any two of the rapidity variables; 

2. /* has the quasi-periodicity property ()5.2|1 : 

3. /* is an entire functions of all its variables; 

reproduces the linear space of fc-particle finite-temperature form factors ff ^{9i, . . . ,9k; L) 
of non- interacting operators \I'. This constitutes a very simple Riemann-Hilbert problem. 
Of course, the solution to this Riemann-Hilbert problem also reproduces the set of zero- 
temperature form factors; in order to disentangle the finite-temperature form factors, one needs 
the additional information about the mixing matrix. 



6.2 Riemann-Hilbert problem associated to right-twist operators 

As recalled in Appendix O there are two types of operators associated to each twist field: those 
with branch cut on their right ("right-twist operators"), and those with branch cut on their left 



21 



("left-twist operators"). As explained in Section |21 with an appropriate choice of the direction 
of the branch cuts, one can obtain two-point functions of twist fields in the NS sector and in 
the R sector. In this subsection, we will consider right-twist operators only. 

One can expect a description for the finite-temperature form factors of twist fields (which are 
"interacting" fields) in the same spirit as the one in the previous sub-section for non- interacting 
fields. 

Consider the function 

f{6i, . . . , Ok, L) = /+,l,+ (^i, . . .,9k;L) 

where is the operator with branch cut on its right representing a twist field: this can be 
the order field (T+ or the disorder field or any of their conformal descendants (that is, fields 
which reproduce conformal descendants in the massless limit). Conformal descendants include 
space derivatives, as well as other fields related to action of higher conformal Virasoro modes 
on twist fields. A way of describing such descendants is by taking the limit x — >■ of the finite 
part of the OPE \E'(x)(j+(0) or \E'(a;)yU_(.(0), where is any bosonic non-interacting field. 

The function / solves the following Riemann-Hilbert problem: 



1. Statistics of free particles: / acquires a sign under exchange of any two of the rapidity 
variables; 

2. Quasi-periodicity: 

f{eu...,ej + 2t7r,...,e,;L) = -f{eu...,e,,...,ek;L), j = i,...,k; 



3. Analytic structure: / is analytic as function of all of its variables 6j, j = 1, . . . ,k every- 
where on the complex plane except at simple poles. In the region Im(^j) G [—27?,^], j = 
1, . . . ,k, its analytic structure is specified as follows: 

(a) Thermal poles and zeroes: f{Oi, . . . ,9^; L) has poles at 



ITT 

%=a„-y, nGZ, j = l,...,k 



and has zeroes at 



ivr 1 
% = On - y, e Z+ -, j = l,...,k, 

(b) Kinematical poles: f{9i, . . . ,9k; L) has poles, as a function of 9k, at 9j ± in, j 
1, . . . , k — 1 with residues given by 

f.n n n I {-l)'~n + e-''=^^ f{9,,...,9„...,9k-i;L) 

f{9l, ...,9k,L)r~u ± -— — _ . 

TT 1 - e 9k - 9j =F m 



In order to have other finite-temperature form factors than those with all positive charges, one 
more relation needs to be used. We have: 
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4. Crossing symmetry: 

f?,:..,e„...,JOu ...,%+ 4; L) = . . . , ^„ . . . , 4; ^) . 

The same relation is also valid for non-interacting fields; recall ()5.1|1 . The name "crossing 
symmetry" is inspired by the zero-temperature case. To make it more obvious, define the 
functions 

f{9[, ...,9[\9,,..., 4; L) = ^_i{9[, 9l\O^{0)\9,, 9^)^,,,^^ . (6.1) 

These are in fact distributions, and can be decomposed in terms supported at separated rapidi- 
ties 9'i 7^ 9j, and terms supported at colliding rapidities, 9'^ = 9j for some i and j. We will 
denote the former by f'^Pi9[, . . . ,9l\9i, . . . ,9k] L), andthe latter by . . . ,9'i\9i, . . . ,9k] L). 
Under integration over rapidity variables, the former gives principal value integrals. Recalling 
the property ()3.14j) . we have 

. . . , 9'i\9i, ...,9k]L) = /++^+ (6'i, . . . ,9k,9'i, . . . , 9[] L) 

for {9[ 9j W i E {1, ...,/}, j G {!,..., k}), where on the right-hand side, there are k pos- 
itive charges (+), and / negative charges (— ). Analytically extending from its support the 
distribution /'^^^- to a function of complex rapidity variables, crossing symmetry can then be 
written 

rPi9[, 9[\9^, ...,9k + i7i]L)= in-{9[, 9[, 9k\9i, 9k-i] L) , 
n-{9[, ...,9[ + i'n\9,,...,9k]L)= in-{9[, 9[_,\9,, ...,9k, 9[] L) , 

which justifies its name. 

It is worth mentioning that the distributive terms corresponding to colliding rapidities satisfy 
a set of recursive equations: 

5. Colliding part of matrix elements: 

f''°'''"{9[, . . ., 9'i\9i, . . . ,9k]L) = 

E m - d,)f{9[, ...,9[,..., 9[\9„ ...,§„..., 9k] L) . 

^ . 1 — e J 

Note that the colliding part vanishes in the limit of zero temperature, L —>■ oo. Finally, it is 
instructive to re- write the distribution f {9[, . . . , 9'i\9i, . . . , 9k] L) as an analytical function with 
slightly shifted rapidities, plus a distribution, using the relation 

^ i7r6{9) + F(]-] (6.2) 



9-i0+ ' ' -\9 



where P means that we must take the principal value integral under integration. Defining the 
disconnected part f disconn. ^g/^^ . . . ,9[\9i, . . . ,9k] L) of the matrix element ()6.H) as 

f{9[, ...,9[\9,,..., 9k] L) = rv-{9[-iQ+, 9[-iQ+\9,, ...,9k] L)+f~{9[, 9l\9,, ...,9k]L) 

(6.3) 
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where again we analytically extend from its support the distribution f^'^P- to a function of 
complex rapidity variables, we find that the disconnected part satisfies the recursion relations 

fdisconn.fnl 0'\Q Q ■ T\ 

/ yPl', ■ ■ ■ , ■ ■ ■ ,Ok,L) — 

I k 

E (1 + '''''''' ) '^(^^ - ■■■A---,m,---Jv---, Ok, L) . 

i=i j=i 

Note that the factor (1 + e'^^'j ) ^(^i - Oj) 

appearing inside the double sum is just the overlap 
+{0i\0j)^, so that the equation above can be naturally represented as a sum of disconnected 
diagrams. 

We will derive below all of the three points in the Riemann-Hilbert problem above, as well as 
crossing symmetry and the recursive formula for the colliding part of matrix elements. 

6.3 Mixing for twist fields 

It is natural to suppose that, like in the case of non- interacting fields, some "mixing" occurs 
between a twist field O and its ascendants - fields of lower dimension and with the same locality 
index and statistics as those of (9 - in calculating the finite-temperature form factors. We 
describe here a conjecture for the way mixing should occur. We will use right-twist operators for 
this description, but we expect that the same mixing should occur amongst left-twist operators. 

First note that the fact that finite-temperature form factors have exponential asymptotic be- 
havior, when all rapidity variables are sent to positive or negative infinity simultaneously, can 
be seen as following from the requirement of convergence of the form factor expansion in 
the region < r < L. Consider a solution /^^ to the Riemann-Hilbert problem above (Points 
1, 2 and 3), with the property that 

+ /?,..., ^fc + /3;L)~(vac|a+(0)|^i,...,^,)e^^ as (3 ±oo 

for some twist operators (9+ of spins s. It will be clear in the next section that such solutions 

exist. Consider a mixing matrix for twist fields whose elements Mq^^ttiL) are parametrized by 

two operators (9+ and This mixing matrix is characteristic of the theory, and its elements 
are non-zero only for operators 0+ that are ascendants of O^. The finite-temperature form 
factors of (9+ are the following linear combinations: 

f{9,, ...,9k;L) = Y, {rnL)f^^ {6^,..., 4; L) . (6.4) 

We do not yet have a full derivation of this, neither a description of the mixing matrix, but it 
seems a natural assumption. 
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6.4 Derivation of the Riemann-Hilbert problem associated to right- 
twist operators 



Point 1. 

The first point of tlie Riemann-Hilbert problem of sub-section 16.21 is a direct consequence of 
tfie definition of tlie finite-temperature form factors and of tlie canonical anti-commutation 
relations of the free fermionic modes. 

Point 3a. 

We now derive the position of the residues and zeroes stated in Point 3a. We first concentrate 
on one-particle finite-temperature form factors. Consider the two-point functions 

g{x,T)= ^(vac|O^(0,0)^^(x,r)|vac)^ . (6.5) 

and 

g{x,T)= ^(vac|V^^(x,r)C^(0,0)|vac)^ . (6.6) 
Their finite-temperature form factor expansions are, from ()5.3j) . 



and 



7,(r t) -- ! dOp'^''^ I l^li^lH JPex-Eer f+^jO^L) ip„^+E„ 



(6.8) 

The one-particle form factors are essentially fixed, up to normalization, by the requirements: 

• convergence of the form factor expansion of g{x, r) ()6.7p in the region — L < r < 0; 

• quasi-periodicity of the analytic continuation in r: g{x,T + L) = —g{x,T) if x < and 
g{x, T + L) = g{x, r) if x > 0. 

For g{x,T), the first requirement is replaced by a convergence in < r < L, and the second 
requirement holds unchanged; they also essentially fix the one-particle form factors. The first 
requirement implies that the form factors f^^{9;L) must go as oc e^^ for some number p as 
9 —>■ oo, and as oc e'^^ for some number q as 6 ^ —oo. The second specifies the positions of 
poles and zeros. 

The second requirement will be satisfied if we can deform the contour in ^-space and take poles 
at appropriate values. With the generalization to multi-particle finite-temperature form factors 
in mind, we will consider these deformations in g{x, r) for a; < only and in g{x, r) for a; > 
only. For x < 0, in (j6.7|) . we can shift the 6'-contour upwards for the term containing e'^^"^ and 
downwards for the term containing e^^"^. Consider the functions 
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The anti-periodicity condition at x < will be satisfied if: the function g+{6) has poles at 
9 = a„ + i7r/2 (recall p.9|) ) for n G Z + the function g-{0) has poles at 6' = a„ — z7r/2 for 
n G Z + i; and e'^'/^g+iO' + i7r/2) + e-'^/'^g_{e' -i7i/2) = for all 9' real (except at the positions 
of the poles). Then, by deforming the contours, we get the sum of 27ri times the residues of the 
poles of g+{0) (equivalently, — 27ri times the residues of the poles of g-{6)) and the resulting 
function is anti-periodic in r. For x > 0, consider (j6.8j) . Define 

~a (B) - ^+'(^'-^) ~n (0) - ^-'(^'-^) 

- ' 9-iO) - Y^^e ■ 

With similar arguments for the function g{x,T) to be periodic: the function g+{6) has poles 
at 6' = a„ — i'n'/2 for n G Z; the function g-{6) has poles at ^ = a„ + z7r/2 for n G Z; and 
e~^^^'^g^{9' — m/2) + e^'^/^g^{9' + m/2) = for all 9' real (except at the positions of the poles). 
There must be no other poles in Im(^^) G [— 7r/2,7r/2] for the four functions g±{9) and g±{9) 
(note that g±{9)/g±{9) are entire functions of 9). Below we will argue that there are no other 
poles in the wider range Im(^) G [— vr, tt]. Assuming further no other types of singularities than 
simple poles, we conclude that g±{9) and g±{9) must have no other singularities than those 
mentioned above in the region Im(^^) G [— 7r,7r]. 

Recall, from the general relation (j3.14p . that /_^+(6';L) = (^f^^{9;L)^ . For the function 
fj^^{9; L), this gives the following conditions: 



fj^^{9] L) has poles at6' = a„ — nGZ and has zeroes at^ = Q;n — y, nGZ 



2' 



fj^^{9] L) does not have poles for \m.{9) G [— tt, tt] except for those mentioned above; 
Re (eT^/+ + (^' - 171/2- L)) = , Re (e'ff^+{9' + z7r/2; L)) = for 9' G 3?. 



The first and second points are the one-particle case of Point 3a. We will verify, in explicit 
calculations of finite-temperature form factors below, that the last point is in fact a consequence 
of the full Riemann-Hilbert problem of sub-section (j6.2p along with (j3.14p . 

The generalization to multi-particle finite-temperature form factors ff^+ 6^.(^15 ■ ■ ■ i^^k'iL) goes 
along the same lines. We need to calculate the multi-point function with k insertions of the 
fermion field ip at k different points. We need to consider two cases: one with these insertions 
on the right of and one with these insertions on its left. We only have to assume that there 
are no poles in |Im(^j — 9j)\ < tt/2 if the two rapidities 9i and 9j are associated to the same 
charge, and in |Im(6'j — 9j)\ < n, except possibly at 9i = 9j, if the two rapidities are associated 
to opposite charges. The poles at 9i = 9j (below, we will show that they are present) are taken 
care of by the prescription stated at the end of Section |21 This prescription does not affect our 
arguments, since we were careful to shift rapidity contours in directions making the two-point 
function convergent at r = with positive difference between the positions of the first and of 
the second operator. Then, repeating the arguments, we find Point 3a. 

Points 2 and 4. 
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Point 2 is evidently a consequence of Point 4. Hence we will show the latter. 

We concentrate again, first, on the one-particle form factors. Consider the function g{x, r) (jfj.Sj) 
and its form factor expansion ()(i.7|l . Consider x < 0, and shift the contour in the first term as 
6^6 + 11^/2 and in the second term as 6 ^ 9 — m jl. This gives, as explained above, a sum 
of residues. The same sum of residues can be obtained, up to an overall minus sign, by taking 
the ^-contour slightly above the line of imaginary part 7r/2 in the first term, and slightly below 
the line of imaginary part — 7r/2 in the second term. This gives again a representation of the 
two-point function valid in the region —L < r < 0, x < 0. From this representation, we can 
shift the ^-contours all the way to the line of imaginary part vr in the first term, and to the 
line of imaginary part — vr in the second term, taking residues of poles, if any. We obtain a 
representation of exactly the same form as ()6.7|) . valid in the same regions of x and r. Since a 
representation of this form is unique, each coefficient of the exponentials should be the same; 
in particular, there should be no poles in the region 7r/2 < lm(0) < tt for the first term and 
— vr < Im(6') < — 7r/2 for the second. The same argument can be applied to the representation 
fl6.8|) of the function g{x, r) in the region a; > 0. The fact that the representation is unique, 
taking into account the prefactor e^^^ in ()6.7|) and ()6.8p and the minus sign occurring when 
going beyond the line of poles at imaginary parts ±7r/2, gives 

which is the one-particle case of Point 4. 

The generalization to multi-particle form factors goes along the same lines. The prescription 
stated at the end of Section El for the case of poles appearing at colliding rapidity variables, 
stays invariant under shifting by ±z7r the rapidity integration lines in the argument above. 

Points 3b and 5. 

We now prove Points 3b and 5. We first concentrate on the two-particle case. Consider the 
distribution 

h{Q,\Q,) ^ {{a{e,)O^i0)a^i9,))), = ^ ^ig^l^^;' . (6.9) 



"2 



where on the right-hand side we use the definition ()6.1|) . Expanding the modes in fermion fields 
using ()4.2j) . we have 



4vr J 



X 



e^((V^(x2)O+(0)^(xi)))z.-2e^((V^(x2)O+(0)V^(xi)))i 

+^e-^((V5(x2)O+(0)^(xl)))z. + e-^((7/>(x2)O+(0)V5(xl)))i 

Consider the first term inside the bracket. We want the main behavior of the expression around 
62 = 61. It is obtained by taking X2 and Xi very large, both with the same sign. In this limit, 
the correlation functions factorize, so the main contribution of the first term is contained into 



^ / dx2rfxie-*P^"^+*P^"^e^sign(-xi)((^(x2)V^(xi)))f"^"^\(O+(0)))L (6.10) 



27 



where the correlation function of fermion fields is taken with periodic conditions around the 
cylinder if Xi > 0, and anti-periodic conditions if Xi < 0: 

((■ ■ ■ trace with periodic (+) / anti-periodic (— ) conditions on the fermion fields . 

The factor sign(— xi) arises because we must put the operator ip{x2) at a slightly positive 
Euclidean time T2 = O"*" and ijj{xi) at a slightly negative Euclidean time ri = 0^ in the initial 
correlation function, and we must take ip{xi) through the branch cut produced by O+{0), if 
xi > 0, before we can factorize. Changing variables, this is 

^ J dx2dxi e-'P'^'+'^'-P'>' e^sign(-xi)((V'(x2)V(0)))?"^"^\(e'+(0)))i . (6.11) 

The integral over Xi is a sum of two distributions, one supported on P2—P1 7^ 0, one supported 
at P2 — Pi = 0. They can be evaluated by using the distributional equations 



/ 



dxe^^'^sign(x) = 2ipQ^ , J dx e'P'^ ^ 27rS{p) (6.12) 



where P means principal value. The part supported on p2 — Pi 7^ is 

p (-^] I dx2 e-^^-^ [mx2)m))t + mx2)m))-L\ {{o^m)L (6.13) 



47r ~ \p2 - pi 
whereas the part supported at p2 — Pi = is 



^ sip2 -Pi) j dx2e-^^^^-e^ [mx2)m))i - mx2)m))i] {{o+m)L . (6.14) 

Putting all terms together, the former gives 

h'-'-me,) - (^^) j dx2e-^^^^^ [moM))t + {{m\x2m] ((o+(o)))l 

(6.15) 

and the latter gives 

/.-"•(^2i^i) = jS{e2 - ^1)-^ J dx2e-^- [{{m\^2m - ((^(^1^2)))^] {{o^m)L 

(6.16) 

where 

-^iOi\x2) = e^i7/'(a;2)^(0) - #(3:2)^(0) + ^^(3^2)^(0) + 6-^^^(0:2)^(0) . 
In order to evaluate the integral over X2, consider the traces 



((a(^2)at(^0))l = 



S{92 - 9-,] 



1 + e-^^oi 

{{a{e2)a\e,)))t = 



A derivation similar to the one above but now applied to these objects gives 



y<5(p2 -Pi) / dx2e-^-- ((*(^ik2)))! = [^^^JeI ■ (6-17) 
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Hence, 



27r-Ve2-^i 



+ 



1 



1 



(l-e-^^''i)(l + e-^^«i; 



(((^+(0)))i 



from which we have 



r^\e2\ei]L) ~ - p 



(((^+(0)))i 



(6.18) 



(6.19) 



Combined with crossing symmetry, this proves Point 3b for the case /c = 2, j = 1. On the other 
hand, the delta-function part is given by 



/^-"■(e^l^l) = ^5(^2 -^l) 



1 



1 + e-^^«i 1 - e-^^»i 



(l_e^^«i)(l + e-^^«i) 
which shows Point 5 for l = k = l,i = j = l. 



(6.20) 



A similar argument holds for k > 2, with extra minus signs coming from the odd statistics of 
the fermion fields, of their modes, and of the twist operator if it has non-zero form factors with 
odd particle numbers only. 



6.5 Riemann-Hilbert problem associated to left-twist operators 

Similar arguments can be used to analyze finite-temperature form factors of twist fields with 
branch cut on their left. We only state here the corresponding Riemann-Hilbert problem. 

Consider the function 

f{9i, . . . ,9k;L) = /+,7.,+ (^i, . . .,9k;L) 
where 0_ is the operator with branch cut on its left representing a twist field. 

The function / solves the following Riemann-Hilbert problem: 

1. Statistics of free particles: / acquires a sign under exchange of any two of the rapidity 
variables; 

2. Quasi-periodicity: 

m,...,ej + 2in,...,ek;L)^-f{0,,...,0j,...,ek;L), j^l,...,k; 

3. Analytic structure: / is analytic as function of all of its variables 6j, j = 1, . . . , k every- 
where on the complex plane except at simple poles. In the region Im(6'j) e [— ztt, ztt], j = 
1, . . . ,k, its analytic structure is specified as follows: 
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(a) Thermal poles and zeroes: /(^i, . . . , (9/^.; L) has poles at 



6'j = a„ + y , n e Z , j = l,...,k 



and has zeroes at 

ITT 1 

= «n + ^ e ^ + 2' J = l'---'^' 

(b) Kinematical poles: f{Oi, . . . ,6k', L) has poles, as a function of 6k, at ± irr, j 
1, . . . , k — 1 with residues given by 

f.n n n ^ {-l)'~^l + e-'^'>^ m,...,6„...,6k-i;L) 



Again, in order to have other finite-temperature form factors than those with all positive 
charges, one more relation needs be used. We have: 



4. Crossing symmetry: 



Moreover, matrix elements 

f{6[, . . ., 6'i\6i, . . .,6k;L) = +^...^{6[, . . . , 6[\0'i{Qi)\6i, . . ., ^fc)^_ _^ 

can again be decomposed in terms supported at separated rapidities 6[ ^ 6j, (which give 
principal value integrals under integration), and terms supported at colliding rapidities, 6'- = 6j 
for some i and j, denoted respectively by f^^^'{6[, . . . ,6[\6i, . . . ,6k;L) and f {6[, . . . ,6[\6i, . . . ,6k 
RecaUing the property (jSHH), we have 

f^''^'{6[, . . . , 6'i\6i, ...,6k;L) = /+~ + „_ . . . ,6k,6'i, . . . , 6[; L) 

for {6[ ^ 6 j 'i i E {1, . . . ,1} , i E {1, . . . , k}), where on the right-hand side, there are k positive 
charges (+), and / negative charges (— ). The distributive terms corresponding to colliding 
rapidities satisfy the same set of recursive equations as in the case of right-twist operators: 



5. Colliding part of matrix elements: 
/™"-(6''^, . . . , 6[\6i, ...,6k;L) = 



1 I — LE, 

1 + e 



E Y.^-ir'~'~' \ LEe ^(^^ - ...A---.0'^6,,...,6„..., 6k. L) . 

■~~r ~T 1 — e J 



i=\ j=i 
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Finally, we can again re-write the distribution /(6''j^, . . . , ^[l^i, . . . , ^fc; L) as an analytical function 
with slightly shifted rapidities, plus a distribution, using this time the relation 

^ = -«W + p(i). (6.21) 

Defining the disconnected part f^i-sco^n.^Q/^^ . . . ,6'i\6i, . . . ,6k] L) of the matrix element ()6.1|) as 

fie[, . . . , e[\e,, ...,9^,1) = r^ie'.+to^, . . . , e[+to^\e,, ...,6^; L)+f~-{e[, . . . , e[\er, ...,6^,1) 

(6.22) 

where again we analytically extend from its support the distribution to a function of 
complex rapidity variables, we find that the disconnected part satisfies the recursion relations 

i-disconn.fnf nf\n O ■ r\ 

I k 

J2 (1 + e-^^*. ) 6{ei - e,)f{e[, ei...,e[\e,,..., §„..., e^-, l) . 

i=i j=i 



7 Finite-temperature form factors of twist fields 

For the order and disorder fields, a and (again, see Appendix O for the definition of the 
associated operators a± and fi±), the solutions to the Riemann-Hilbert problem of sub-section 
16.21 are completely fixed (up to a normalization) by the asymptotic behavior ~ 0(1) at — > 
±00, since they are primary fields of spin 0. Note that the method of computing one-particle 
finite-temperature form factors by solving the Riemann-Hilbert problem with this asymptotic 
behavior is very similar to the method used by Fonseca and Zamolodchikov j27j for calculating 
form factors on the circle. 

For the one-particle finite-temperature form factor of the disorder operator with a branch cut 
on its right, the solution is 

f^^e;L) = e^fC{L) exp 

for some real constant C{L). This is in agreement with the Hermiticity of which gives 
(/^+(^; L))* = f^+{e; L) for 6 real. Using 

sinh(^ - {9' ± 20+)) " ^'""^^^ - + P (gj^j^^^ _ ■ 
this can also be written 

/ 1 _|_ p-LEg 

(7.2) 



de' 



oo=F^0^ 



27Ti smh{e - 6'] 



In 



l + e" 



-LEal 



1 - e 



(7.1) 



J_^ 2-ni- V sinh(e -9')) ^ U - e-^^o' ) 
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That this is a solution can be checked by verifying the asymptotic behavior f'^{6]L) ~ 
e^^C{L) as \0\ oo, and by verifying that the functions f'^{9]L) have poles and zeros 
at the proper positions. Positions of poles and zeros are the values of 9 such that when analyt- 
ically continued from real values, a pole at sinh(^ — 9') = Q in the integrand of ()7.1|) and one of 
the logarithmic branch points pinch the 9' contour of integration. The fact that these positions 
correspond to poles and zeros can be deduced most easily from the functional relation 

1 I -LEe 

f^^9;L)ft^9 ±trr;L) = ±tC{LY^^^ . (7.3) 

Note that this implies the quasi-periodicity property 

f^^9 + 2zn-L) = -f^^9;L) . (7.4) 

It is also easy to see that the crossing symmetry relation is satisfied. Also, since C{L) is real, one 

can check the validity of the relation Re ( e~ f_^_^ {9' — iTT/2] L) ] = , Re ( f_^_'^{9' + i'K/2; L) 

for 9' real; this relation was seen as a consequence of general principles in the proof of Point 
3a in sub-section 16. 4j it is now seen as a consequence of the Riemann-Hilbert problem along 
with Hermiticity. 

For the operator /i_ with a branch cut on its left, one can check similarly that the function 

r±- {9- L) = {9 - L) = -if^^ {9- L) (7.5) 
solves the Riemann-Hilbert problem of sub-section 16.51 Explicitly, 



f^-{9-L) = -ie^fc{L) exp 



d9' 1 , fl + e-^^" 
± / ^-7z TV In 



oo±jo+ 27ri sinh(6' — 9') \1 — e 



-LE, 



(7.6) 



In particular, we observe that {f±~ {9; L))* = —f^~ {9; L), which is in agreement with the anti- 
Hermiticity of the operator fi^ (see Appendix E)). Note that we chose the same constant C{L) 
as a normalization for both fl^~ and /j^^. This is not a consequence of the Riemann-Hilbert 
problem, but can be checked by explicitly calculating the normalization. The normalization is 
calculated in Appendix |Dl and is given by 

C(L) = (7.7) 



where the average {{<j))l was calculated in ^29^ (the average at zero-temperature (that is, 
L — s> oo) can be found in j3B]) and is given by 



1 J_ _ 1 3 

m8 2i2e sA^exp 



mLf f d9id92 sinh0isinh^2 

' ' In 



2 J J _^ (271")^ sinh(mL cosh 6'i) sinh(mL cosh ^; 
where A is Glaisher's constant. 



coth — — ^ 



Multi-particle finite-temperature form factors can be easily constructed from the well-known 
zero-temperature form factors (first calculated in [H]), by adjoining "leg factors", which are just 
normalized one-particle finite-temperature form factors: 
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where 0+ is 0"+ if k is even, and /i+ if k is odd. The symbol [k/2] equals the greatest integer 
smaller than or equal to k/2. This satisfies the condition on thermal poles and zeroes simply 
from the properties of the leg factors, and it can be verified that this satisfies the quasi- 
periodicity condition and the kinematical pole condition. Point 2 and Point 3b of sub-section 
l(i.2| respectively. Using crossing symmetry. Point 4, it is a simple matter to obtain the formula 
for other values of the charges: 

Similarly, we have 

/,t„Jft. .... ft; L) = a« ((.)). (n ^^5y^) ^ (^) ) 

where (9_ is (T_ if k is even, and yU_ if k is odd. 

It is easy to check, using (|1.8p . that the formulas above for finite-temperature form factors 
reproduce the known form factors on the circle [2^1 12ni HZ] • 

Also, it is a simple matter to obtain a Fredholm determinant representation for the two-point 
function of twist fields in the R sector. This is derived in Appendix |E| 

The zero-temperature form factors of descendant twist fields in the Majorana theory were 
analyzed in |17j, as solutions to the form factor equations. It was found that the set of form 
factors of descendant fields of a given spin s is described by multiplying form factors of primary 
twist fields by symmetric polynomials in e^^ , j = 1, . . . , (where 6j are the rapidities of the form 
factor) homogeneous of degree s. Similarly, we expect that finite-temperature form factors of 
descendant twist fields are obtained by multiplying those of primary twist fields by symmetric 
polynomials. It is easy to check that the requirements of the Riemann-Hilbert problems of 
Sections 16.21 and 16.51 are still satisfied after such an operation. However, as explained in the 
previous section, since rotation is no longer a symmetry on the cylinder, we cannot identify 
the spin of a descendant with the homogeneous degree of symmetric polynomials. Instead, we 
must rely on the conjecture in sub-section 16.31 in order to describe how to combine appropriate 
symmetric polynomials. Unfortunately, we have not yet calculated the mixing matrix for twist 
fields, hence we cannot derive more explicit formulas here for descendants. 

8 Conclusions 

We have derived a method for writing large-distance expansions of finite-temperature corre- 
lation functions. For this purpose, we introduced the space of operators C on which finite- 
temperature quantum averages are vacuum expectation values, and we obtained the expansion 
of two-point functions by inserting a resolution of the identity between the two operators. We 
defined finite-temperature form factors, which are appropriately normalized matrix elements on 
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L of local operators, and we described at length their properties. They are related, by analytical 
continuation in rapidity space, to form factors in the quantization scheme on the circle. Any 
finite-temperature form factor of a non-interacting field can be seen as the zero-temperature 
form factor of a sum of non-interacting fields, including itself and fields of lower dimension. 
The coefficients in this sum are independent of the matrix element evaluated. We described 
this phenomenon by constructing the associated mixing matrix; note that this was overlooked 
in the initial approach 31, . We derived the Riemann-Hilbert problem that defines the set of all 
finite-temperature form factors of twist fields, and calculated finite-temperature form factors 
of order and disorder fields. 

There are still many open questions: 

We have only implicitly described how to obtain large-distance expansion of correlation func- 
tions with more than two operators: one should insert the resolution of the identity between 
every pairs of adjacent operators. The resulting expression involves matrix elements where both 
vectors in L correspond to excited states. For non-interacting fields this does not cause any 
difficulties, but for twist fields, one must be more careful. We evaluated such matrix elements 
of twist fields, as distributions with terms supported at colliding rapidities (which vanish at 
zero temperature) and others supported at separated rapidities. Using these distributions, all 
integrals over rapidity variables are well-defined, and it would be nice to explicitly observe the 
agreement of the resulting expression for multi-point correlation functions with a form factor 
expansion on the circle. In particular, from the results presented here, it should be possible to 
calculate matrix elements of twist fields in the quantization on the circle where both vectors 
correspond to excited states. 

We have described how to obtain form factors on the circle from finite-temperature form factors 
only for excited states in the NS sector, where fermion fields have anti-periodic condition around 
the cylinder. It would be interesting to describe matrix elements in a similar way with excited 
states in the R sector. As we mentioned, this can be done by defining the inner product on L 
as a trace with a twist operator at, say, position x = — oo and having a branch cut on its right. 

We have explained how to obtain the mixing matrix associated to non-interacting fields. In 
particular, we have obtained an operator U that describes the action of this matrix on the 
non-interacting fields seen as vectors in C. The existence of such an operator is well-known 
in conformal field theory: there, the space of operators is isomorphic to the Hilbert space in 
radial quantization (which is isomorphic to the Hilbert space in any quantization on closed 
lines), and the mixing operator is described by an operator on the Hilbert space performing a 
transformation to the cylinder. It would be interesting to relate more explicitly the operator 
U that we wrote with a transformation to the cylinder. 

We obtain finite-temperature form factors of order and disorder fields. It is important to 
note that form factors for the right- and left-operators both specialize to the zero-temperature 
form factors only in the region |Im(^)| < 7r/2 of the complex plane of rapidity variables 6. 
In connection to this, observe that the quasi-periodicity relation does not agree with the one 
satisfied by zero-temperature form factors. At finite temperature, this quasi-periodicity relation 
does not seem to contain any information about the mutual locality of the twist field with 
respect to the fundamental fermion fields. This semi-locality is included in the thermal poles 
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and zeroes and in the particular form of the kinematical residue. It would be interesting to 
generalize this to other semi-locality index, by deriving the Riemann-Hilbert problem in the 
free Dirac theory for scaling twist fields associated to the U{1) symmetry (in this connection, 
see the recent results jlH])- But the thermal poles and zeroes might not be sole consequences 
of semi-locality; they might be interpreted as coming from self-interaction of the field around 
the cylinder as well, hence could be present also for fields with zero semi-locality index in 
interacting theories. 

It could interesting to obtain the non-linear differential equations describing the two-point 
function of twist fields from the Fredholm determinant representation resulting from the finite- 
temperature form factor expansion. This would modify some of the results of j^; in particular, 
it is easy to see that the equation evolving temperature derivatives would not hold due to the 
non-trivial leg factors, and it is not clear if there is an equation replacing it. 

Finally, probably the most important future development from the ideas presented here is the 
definition and calculation of finite-temperature form factors in interacting integrable models, 
and the identification of the associated measure of integration in the resolution of the identity. 
An important new idea presented here is the relation between finite-temperature form factors, 
seen as matrix elements on the Hilbert space of operators, and form factors on the circle; in 
particular, the relation between the measure of integration and the spectrum on the circle. If, 
for fields that are mutually local with respect to fundamental fields, the finite-temperature form 
factors equations are the same as the well-known zero-temperature form factor equations, as was 
assumed in previous approaches, then the only non-trivial element is the measure of integration 
in the resolution of the identity. This measure was indeed subject of controversy [SSI EH 
Our ideas give a guideline for defining this measure of integration. Note, however, that in view 
of the works 1221; rapidity variables might turn out not be the most appropriate ones 
for defining finite-temperature form factors. 
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A Derivation of the relation ( 11.81 ) 



A general argument for formula ()1.8|) can be obtained from the relation between traces and 
expectation values on the circle, (jl.lj) and (jl.3|) . Recall the expression of mode operators in 
terms of fermion operators (j4.2j) . From this, we can write 

= {{a{9l)■■■a{9,)O{0,0)a\e^)■■■a\9k)))L (ri(l + e" 

\j=i 
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where we define the operators 

We are interested in the analytical continuation 6 d + i7T/2 of the matrix element ()A.1|) for all 
rapidity variables, then in taking the limit where 6 ^ an ()1.9p of this analytical continuation. 
The result of this limit is obtained by taking the analytical continuation W±{6) = w±{6 + i7r/2) 
of the operators above, then by taking, in all operators W^{6) on the right of 0(0), only the 
negative part of the integral over x, and in all operators W-{6) on the left of 0(0), only the 
positive part of the integral over x. Indeed, the variable x is a time variable in the quantization 
on the circle, positive x corresponding to positive time, and the result of the limit is obtained 
by looking at the time-ordered part of the correlation function. 

Hence we have 

W+{e) ~ (l + e-'^Po) / (ixe-^''He5+fV^i(0,x) -ie~^-'^7/^L(0,x) 

J —00 

f,-Eex 



(1 + e-'^P") 



dx 



A/cosh(a„) 



X 



e ( e 2 + e 2 



an + e-^"^ I e 2 - e 



(1 + e^'^P") 



E 



1 



2L , \/ coshfon) 
36 



e2+e 2 62— e 



En + -Eg 



In the last step, we used analytical continuation in the exponents in order to perform the 
integral. We can now extract the result of the limit 9 —>■ a„: 

w+{e) ~ 



i 2sinh(a„) 1 + e-'^P» . 
—— — d ' 

2L A/cosh(a;„) En — 



a/2L cosh(a„) ajj 



(A.3) 



On the other hand, 

^Eex 
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dx 
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e ( g 2 _ g 2 ) a„ + e^"^ ( e 2 + e 2 ] 
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(1 + e-'^P') 



E 



which gives the result of the limit 6' — > a„ to be 

W_{9) ~ V2Lcosh(a„) a. 

Hence we recover ()1.8|) . 
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B Proof of the equivalence between the finite-temperature 
form factor expansion and the form factor expansion 
on the circle 



By re-arranging the rapidity variables (along with their associated charges) in p.fjp . it is possible 
to bring the finite-temperature form factor expansion of two-point functions in the form 



X 



00 K 

^oi^oJ k\{K- k)\ n ■=! (1 + e-™^-^*^^.) nf=.+i (1 + e™^^°^^^^ 

X /+^^..^+ _(6'i, . . .,9k, 9k+i, . . . , 9k] L)f^^_^__^_{9K, • • • , 9k+i, 9k,..., 9i; L) . 

(B.l) 

If the first field Oi is a twist field, we will consider it to have a branch cut on its right (positive 
X direction); on the other hand, if the second field O2 is a twist field, we will consider it to 
have a branch cut on its left (negative x direction). This ensures that when we obtain the 
form factor expansion on the circle, the intermediate states are in the sector with anti-periodic 



37 



conditions on the fermion fields (whereas the vacuum vectors may be in the sector with periodic 
conditions). Moreover, we will assume a; > 0; this insures that the operators are time-ordered 
in the quantization on the circle. 

Consider shifting the contours of integration associated to the rapidity variables 6i,...,6k 
towards the positive imaginary direction, by an amount in. The only poles that contribute are 
those from the factors 11^=1 (l + ^-^LcoshOj^ ^^^e denominator. 

Consider the terms arising from taking poles, with fixed K. They are given by 
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( %Tl ZTT 

X f_^l^^^j^_^^^ _ \ Ok-, • • • , On+ii ttnjv + ^5 • • • 5 + 

where each element of the set {rij, j = 1, . . . , A^} is a number in Z + i. On the third line, there 
are positive charges and K — N negative charges; on the fourth line, there are K — N positive 
charges and negative charges. On the second line, the factor (— 1)'^"^ comes from shifting 
k — N rapidity variables by ivr, and taking the imaginary factors from the crossing relations, 
()5.1|) and Point 4 in sub-section 16.21 The sum over k can be done, and vanishes whenever 
N ^ K. Hence we are left with 

.Ef=i(-"^^cosha„^^-n,2Eir) ^ 



n 



fO ( ^'^ t\ ro-, ( in ^ 

X /+,!..,+ ( "ni + y , • • • , + y ! ^ ) .. - 1 "n/f + • • • ' "ni + y ; 

When summed over K and over {rij}, this reproduces the form factor expansion on the circle 
as in ()1.4|) if we use ()1.8|) . up to a phase factor in accordance with ()1.3j) . 



C Operators on H corresponding to twist fields 

Let us now turn to the description of the operators associated to twist fields: the order field a 
and disorder field /i. 

The order and disorder fields can be essentially defined by their OPE's with the fermion fields. In 
order to keep the proper OPE's, it is convenient to remember that with the conformal normaliza- 
tion 2;') ~ (2;— the leading term of the OPE's are '?/'(2;)a(0) ~ \/^/^(0)) '?/'(-s)/^(0) ~ 
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■2^o'{0). Our normalization of the fermions, however, is given by ()2.9p . Hence we have 

^(x, r)a(O) ~ , / ^ . /i(O) , ^(x, r)^(O) ~ , / ^ . a(0) (C.l) 

and 

^(x, r)a(O) ~ - , / . ^ ^(0) , ^{x, r)/i(0) ~ , / . ^ a(0) . (C.2) 

By convention, the disorder field has nonzero odd-particle form factors only, and the order 
field has nonzero even-particle form factors only. In particular, a has a nonzero real vacuum 
expectation value. 

The finite-temperature correlation functions of fermion fields with insertion of these spin fields 
are functions on coverings of the cylinder. In order to define operators on the Hilbert space Ti 
corresponding to the twist fields a and fi, we must choose Riemann sheets. In fact, since the 
quantization is on the line, the cylindrical geometry is implemented by tracing over Ti, so that 
we may define the operators corresponding to twist fields by specifying the Riemann sheets on 
the plane (that is, it is sufficient to define them at zero temperature). It will be convenient 
to define two operators on 7i for each twist field, corresponding to different choices of branch 
cuts. Inside vacuum expectation values in 7i, the choice of branch cut does not matter, but 
inside traces, it does. 

Consider the product of fields {ip{x, r)yu(0))+ defining (inside zero-temperature correlation func- 
tions) a function on the plane satisfying the free massive equation of motion everywhere except 
at the branch cut r = 0, x > (and except at the position of other local fields, if any). We will 
define the operator /i+ by the fact that under the mapping of matrix elements of operators to 
correlation functions, 

V'(x,r)/i+(0) ^ (^(x,r)MO))+ (r > 0) 

/x+(0)^(x,r) ^ -(V^(x,r)/i(0))+ (r < 0) (C.3) 

and 

^(x,r)^+(0) ^ [(V^(x,r')^(0))J (r < 0) 

/x+(0)^(x,r) ^ -[(^(x,r')/i(0))+l (r > 0) . (C.4) 

That is, in the second set of maps, the product of operators gives, inside vacuum expectation 
values, correlation functions that are continued through the branch cut of (■?/'(x, r')/i(0))+ if 
X > 0. 

Similarly, consider the product of fields ('?/'(x, r)/i(0))_ defining a function on the plane that 
satisfies the free massive equation of motion everywhere except at the branch cut r = 0, x < 
(and except at the position of other local fields, if any); this function coincides with the function 
defined by ('?/'(x, r)/i(0))+ when r > 0. We will define the operator /i_ by the fact that under 
the mapping of matrix elements of operators to correlation functions, 

V'(x,r)/i_(0) ^ (^(x,r)MO))- (r > 0) 

/i_(0)^(x,r) ^ -(^(x,r)/i(0))_ (r < 0) (C.5) 
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and 



?/'(x,r)/i_(0) 
^_(0)V^(x,r) 



(^(x,r')/x(0)) 
(^(x,r')M(0))_ 



r'=0+- 



r'=0-^T 



(r<0) 
(r>0) 



(C.6) 



We make similar definitions for the order field, without the minus sign that represents the odd 
statistics of fi± with fermion operators. We have: 



V(x,r)a+(0) ^ (V^(x,r)a(0))+ (r > 0) 
a+(0)^(a;,r) ^ (^(x, r)(T(0))+ (r < 0) 



(C.7) 



and 



V^(x,r)a+(0) 
a+(0)^(x,r) 



(V^(x,r')a(0))_ 
(^(x,r')a(0))_ 



=0+- 



(r<0) 
(r>0) 



(Ci 



as well as 



?/'(x,r)(T_(0) (^/'(x,r)a(0))_ (r > 0) 

(T_(0)V^(x,r) (^(x,r)cr(0))_ (r < 0) 



(C.9) 



and 



^(x,r)a_(0) 
a_(0)V^(x,r) 



(7A(x,r')a(0))_ 
(V^(x,r')a(0)) 



=0+- 



(r<0) 
(r>0) 



(C.IO) 



With these definitions and with the OPE's ()C.1|) and ()C.2|) . it is possible to check that the 
fields a± are Hermitian, that /i+ is Hermitian and that yU_ is anti-Hermitian on Ti: 



(C.ll^ 



Note that relations ()C.3I - IC.10|) also hold for operators corresponding to descendants of the 
order and disorder fields. 



D Normalization of the one-particle finite-temperature 
form factors 

The normalization can be obtained by considering the OPE's tp{x,T)fi±{0). Consider first 
/i_(0). We have 

((rz^(x,rK(0)))L ~ , / , . A {a))L (D.l) 
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where T means time ordering (the latest operator being placed on the left) and where the 
square root is taken on its principal branch. We will take r > 0, x > in ()D.1|) : then we 
can repeat the calculation of Appendix ^ Since only the one-particle finite temperature form 
factors contribute, the results of Appendix iBl give, for the left-hand side, 

' mLcoshK)^H J r" + T'^J • 

The leading x + ir ^ 0^ behavior is obtained by looking at the leading |n| — > oo behavior of 
the summand. We have 

r-- [a^ + Y,L^-e-'^C{L). 

From ()5.3|) . the main contribution comes from the region — > cxd only, and we obtain for the 
leading behavior of the left-hand side of ()D.1|) . 



nGZ+i,n>0 



Hence, we find fl7.7p . 

It is instructive to repeat the calculation for /i+(0) and to verify that C{L) is indeed as given 
above. We have 

((r^(a:, r)/i40)))i ~ , , / . (D.4) 

where the square root is on a branch that coincides with the principal branch for r > but 
that has a cut at r = 0, a; > 0. In order to apply the results of Appendix |B1 we need take 
T < 0, X < 0, and we have 

((/i+(0)^(a;,r)))^ ~ t ] , {{a))^ (D.5) 

2a/— vr(x + tT) 

where now the square root is on its principal branch. This has expansion 
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With f1^{6] L) ~ e^i C{L) and using we have, for the leading x + ^ behavior of 

the left-hand side, 



iC{L) 27rn(x + iT) 1 



L ^ ^ \-2(x^iT) 

n€Z+^,n>0 v v / 



which agrees with ()D.5|) if C{L) is given by ()7.7|) . 
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E Fredholm determinant representations for two-point 
functions at finite temperature 



The results of Appendix iBl essentially show that 

(E.l^ 



r)(T_(0,0)))L = L(vac|crL(-r,x)crL(0,0)|vac)L 



r);U_(0,0)))L = L(vac|;UL(-r,x)/iL(0,0)|vac) 



L 



where on the right-hand side, the vacuum expectation values are on the circle in the Ramond 
sector. In a similar fashion, an analysis of the finite-temperature form factor expansion shows 
that 

((a_(x,r)a+(0,0)))L= (K(-x,-r)a_(0,0)))i 
((/x_(x,rK(0,0)))L = -((/i+(-x,-r)/x_(0,0)))L. ^ ' ^ 

Note that the second equation is in agreement with the fermionic statistic of the operators 

Using the finite-temperature form factors ()7.9j) and ()7.10|1 for twist fields, we have then the 
following large distance expansions of two-point functions in the R sector: 

((a+(a;,r)a_(0,0)))i = 

°° r rIO AO e, (imzsinh Si— mr cosh 0, ) ^ / O O \ 

k even 

(E.3) 

and 



((/x+(x,r)/i_(0,0))) 
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y y / -^'WUr{Of,L)f W tanhf^ 

k=0 ti,...,ek=±-^ \.\.j = l\^^^ I j = l l<i<j<k ^ 

k odd 

(E.4) 

Following |M] , Fredholm determinant representations can now easily be obtained from the 
formulas 



detjj- 



k odd 



and 



^^^''^ { Mi + } . . . Uk i<n^fc ^Ui + uj) ' ^^'^^ 

Formula ()E.5|) gives 

(((T+(x,r)a_(0,0)))i = det(l + K) (E.7) 

where K is an integral operator with an additional index structure, defined by its action 
(K/)e(^) = Ee'=±IZodO'K,,AS,0')fAO') and its kernel 

(n/ n\ e{imx sinh 8— niT cosh 8) 

^) 1 + e-^^.cosh. ■ (E-8) 
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Finally, in order to obtain two-point functions of disorder fields, we must consider the linear 
combinations a ± fi. Formula ()E.6|) gives 

{{{a4x, r) + r//i+(x, r))(a_(0, 0) + r7/i„(0, 0))))l = det(l + (E.9) 

with 7] = ± and by definition {J^'^^f)^{u) = du' J^'}i{u,u')f^i{u') where the kernel is 

given by 

1 p^{imx{u—u~^)—mT(u+u~^)) 
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